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Abstract 

In this paper, we revisit the forward, backward and bidirectional Bahl-Cocke-Jelinek-Raviv (BCJR) 
soft-input soft-output (SISO) maximum a posteriori probability (MAP) decoding process of rate-1 
convolutional codes. From this we establish some interesting duality properties between encoding 
and decoding of rate-1 convolutional codes. We observe that the forward and backward BCJR SISO 
MAP decoders can be simply represented by their dual SISO channel encoders using shift registers in 
the complex number field. Similarly, the bidirectional MAP decoding can be implemented by linearly 
combining the outputs of the dual SISO encoders of the respective forward and backward decoders. The 
dual encoder structures for various recursive and non-recursive rate-1 convolutional codes are derived. 
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I. Introduction 

Convolutional codes were first introduced by Elias more than 50 years ago [1]. They have 
been widely used in various modern communications systems, such as space and satellite com- 
munications, cellular mobile, digital video broadcasting etc. Its popularity stems from its simple 
encoder structure, which can be implemented by using shift registers. 

The main complexity associated with systems using convolutional coding is situated in the 
decoder. Decoding essentially consists of finding an optimal path in a trellis based graph. Various 
decoding algorithms have been developed to achieve the optimal decoding performance in the 
most efficient manner. The Viterbi algorithm (VA) has been known as a maximum-likehood (ML) 
decoding method, which minimizes the sequence error rate [2-4]. It exhaustively searches all 
states of the trellis over a fixed length window and finds a most likely information sequence. In the 
standard VA, the decoder produces hard-decision outputs, which are the estimates of transmitted 
binary information symbols. In [5, 8], the VA is modified to deliver not only the most-likely 
binary signal sequence, but also the soft output containing the a posteriori probabilities (APPs) of 
the transmitted binary symbols. The soft-output VA (SOVA) is especially useful when decoding 
concatenated codes, such as turbo codes, as it provides soft input for the next decoding stage 
and thus improved performance. 

There exists another class of non-linear decoding algorithms, called maximum a posteriori 
probability (MAP) decoding. It was first proposed by Bahl, Cocke, Jelinek and Raviv (BCJR) 
in 1974 [6]. It performs symbol by symbol decoding and uses the symbol error rate as the 
optimization criterion. Both the input and output of the decoder are soft information signals. 
Compared to the VA, the soft-input- soft-output (SISO) MAP can provide the optimal symbol- 
by-symbol APP, and thus can fully exploit the full benefits of soft-decision decoding in iterative 
decoding process of concatenated codes. 

The BCJR MAP decoding is a bi-directional decoding process, consisting of a forward and 
a backward recursion process, which dominates the main complexity of a decoder. In each 
direction, the decoder infers the probabilities of current states and information symbols based on 
the probabilities of the previous states in the forward and backward trellis, the received signal, the 
channel state and the a priori probabilities of the transmitted signals. The complexity of forward 
and backward recursion exponentially increases with the constraint length of convolutional codes. 
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In this paper, we revisit the forward, backward and bidirectional BCJR SISO MAP de- 
coding of rate-1 convolutional codes. We observe some duality properties between a SISO 
forward/backward BCJR MAP decoder of a convolutional code and its encoder. The forward and 
backward decoder of a rate-1 convolutional code can actually be represented by its corresponding 
dual encoder using shift registers in the complex number field. This significantly reduces the 
the original exponential computational complexity of MAP forward and backward recursion to 
the linear complexity. Similarly the bidirectional BCJR MAP decoding can be implemented 
by linearly combining the outputs of the dual SISO encoders of the respective forward and 
backward decoders. With logarithm of the soft coded symbol estimate, directly obtained from 
the received signals, as the input to the dual encoder, the dual encoder output produces the 
logarithm of the soft symbol estimates of the binary information symbols. We found that the 
dual encoder structure of a code depends on whether the code is recursive or not. We divide 
the rate-1 convolutional codes into three classes, feedback only convolutional (FBC) code, feed- 
forward only convolutional (FFC) code and general convolutional (GC) code. The dual encoder 
structure is derived for each class of codes. 

The remainder of the paper is organized as follows. In Section II, we first briefly review 
the BCJR forward decoding algorithm and derive the dual encoder structures of BCJR forward 
decoders for three classes of rate-1 convolutional codes. The duality for backward decoding is 
presented in Section III. The representation of bidirectional MAP decoding by using the derived 
dual encoder structures of forward and backward decoding is described in Section IV. Simulation 
results are shown in Section V. Conclusions are drawn in Section VI. 

II. Duality of Encoding and Forward Decoding of Rate-1 Convolutional 

Codes 

In this section, we first revisit the forward BCJR MAP decoding algorithm. We will focus 
on the decoding of a single constituent convolutional code of rate-1. Let b = (61,62, •• ■ ,bx) 
be a binary information symbol sequence to be transmitted, where K is the frame length. Let 
c = (ci, c 2 , . . . , c K ) be the binary codeword of b, generated by the binary code generator 
polynomial g, and x = (xi, X2, ■ ■ ■ , xk) be the modulated symbol sequence of c. For simplicity, 
we consider a BPSK modulation. Let y = (2/1,2/2, • • • ,Vk) denote the received signal sequence 
at the output of channel. 
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Based on the encoder structure, we define three different classes of convolutional codes. Let 
a(x) = x n + a„_ix n_1 + • • • + a\Xi + 1 and q(x) = x n + g n _ 1 x n ~ 1 + • • • + q\Xi + 1, where n 
is called the degree of polynomials a(x) and g(x). We define a convolutional code, generated 
by 9fbc( x ) = l/q{x), as a feedback-only convolutional (FBC) code, a code generated by 
9ffc( x ) — a ( x ) as a feed-forward only convolutional (FFC) code, and a code generated by 
9gc( x ) — a ( x )/g( x )> as a general convolutional (GC) code. We will investigate the forward 
decoding process of these three classes of convolutional codes. 

A. Forward decoding of a FBC code 

In this subsection, we first investigate the forward decoding of an FBC code. To gain better 
insight into the decoding process, let us first look at the following example. 

Example 1: We consider a FBC code with the generator polynomial of gFBc( x ) — x s+ x+ i > 
for which the encoder and trellis diagram are shown in Fig. [IJ In the trellis diagram, the state 
is labeled as Si 62 , where Si, i = 1,2 is the value of the i-th encoder shift register content. 
Each branch in the trellis is labeled as x/y where x and y denote the encoder input and output, 
respectively. 

Let p Ck (l) = p(ck = l\yk), I = 0,1, denote the a posteriori probabilities (APP) of the encoded 
symbol Ck = I, given the received signal y^, where Ck is the transmitted binary coded symbol at 
time k. Let us further denote P c = {(p cl (0),p Cl (l)), ■ ■ •, (p Ch (0) , p Ch (l)) , (p CK (0),p CK (l))}. 
Now let us follow the BCJR MAP forward decoding algorithm to use P c to calculate the 
APPs of binary information symbols bk- Let Pb k {w) = p(bk = w\y) represent the probability of 
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information symbol b k = w, w=0, 1, given the received signals y = {y 1 , ■ ■ ■ , y k , ■ ■ ■ , y K }. It 
can be calculated in the following recursive way 

p bk (w) = p(h = w\y)= ^2 a jfc _i(m / )7fc(w'm) (1) 

(m',m)eU(b(k)=w) 

= a k -i(m')p Ck (ck(Tn',m)) 

(m',m)€U(b{k)=w) 

®k{m) = ^a fc -i(m / )7 fc (m / m) = ^ ak-i(m')p Ck (c k (m' , m)) (2) 

ml m' 

where U(b(k) = w) is the set of trellis branches from the state m! at time k-\ to the state m 
at time k, that are caused by the input binary symbol b(k) = w, and Ck(m',m) represents the 
encoder output of the corresponding trellis branch. 

Let m = 0,1,2,3 represent the states of Si £2 — 00,01,10,11 at time k, respectively. Let 
x c = (x Cl ,--- ,x CK ) and Xb = (^fen - '' i^b K ) denote the soft symbol estimate sequence of 
codeword c and information sequence b, respectively. We assume that and 1 are modulated 
into the symbol 1 and -1. Then the soft symbol estimates x Ck and Xb k , which represent the 
probabilistic average of estimates of symbols x Ck and Xb k given y, can be calculated as 

x Ck = E{x Ck \y k ) = p Ck (0) - p Ck (1) (3) 
%b k = E(x bk \y)=p bk {0)-p bk {l). (4) 

Then by using Eqs. (OQ) and © alternatively in Example 1, we get 

(1) at time k = 0, 

Qt (0) = 1; Qt (l) = 0; a {2) = 0; a (3) = 0; 
p 6o (0) = l;p 6o (l)=0; 

(2) at time k = 1, the received signal is y(l), and the input to the decoder is the APPs of c\, 
given by p Cl (0) and p ci (l), respectively. Then we have 

ai(0)=p Cl (0); ai(l) = 0; ai(2) = p Cl (l); ai(3) = 0; 

p&i(o) =p C i(o); ^(1) =Pci(i); 

and 

Xh = Ph (0) - p bl (1) = Pci (0) - Pd (1) = \x ci ] 

l l 

(3) at time k = 2, the input to the decoder is the APPs of C2, p C2 (0) and p C2 (l). We have 
« 2 (0) =p C2 (0)p Cl (0); a 2 (l) =p C2 (0)p Cl (l); a 2 (2) = p C2 (l)p Cl (0); a 2 (3) = p C2 (l)p Cl (l); 
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VbM =p ca (0)a 1 (0) + p C2 (l)a 1 (2); p b2 (l) = p C2 (l)ai(0) +p C2 (0) ai (2); 
and 

=P6 2 (0) -P6 2 (l) = (Pc 2 (0) -p C2 (l))(p cl (0) -p ci (l)) = ^X C2 Z C J 

1 1 

(4) Similarly, at time 3, we have 

as(0) =p C3 (0)p C2 (0); as(l) =p C3 (0)p C2 (l); a 3 (2) = p C3 (l)p C2 (0); a 3 (3) = p C3 (l)p C2 (l); 
Pb .M =p C3 (0)a 2 (0) +p C3 (l)a 2 (l) +p C3 (l)a 2 (2) + p C3 (0)a 2 (3); 
p 63 (l) =p C3 (l)a 2 (0) +p C3 (0)a 2 (l) +p C3 (0)a 2 (2) + p C3 (l)a 2 (3); 
and 

Xb 3 =Pb-M -Pb 3 (l) = (pca(O) -p C3 (l))(Pc 2 (0) - Pc2 (l))( Pci (0) -p Cl (l)) =^ C3 x C2 x c J 

1 1 

At any time instant k > 2, we can generalize that 

MO) = PcMPc k _M; a k (l) = ^(0)^(1); a k (2) = ^(1)^(0); a fc (3) = ^(1)^(1); 
P6 fc (0) = p Cfc (0)a fc _!(0) +p Cfc (l)a fc _i(l) +p Cfc (l)a fe _ 1 (2) +p Cfc (0)a fc _i(3); 
P6 k (l) = p Cfc (l)a!fc-i(0) +p Cfc (0)« fe _i(l) +p Cfc (0)a fc _i(2) +p Cfc (l)a fc _i(3); 
and 

= Pf. fc (0) -P6 fe (l) = (p Cfe (0) -p Cfc (l))(a fc _i(0) + a fc _i(3) - a fc _i(l) - a fc _i(2)) 

r _ . _ j 

I ry iy /y I 

— 1 » t 'C fe » t 'C fc _l J 'Cfc_ 2 I 
I I 

where 

(a fc _i(0)+a fc _i(3)-a fc _i(l)-a Jfc _i(2)) = p Cfc _ 1 (0)p Cfc _ 2 (0)+p Cfe _ 1 (l)p Cfe _ 2 (l)-p Cfc _ 1 (0)p Cfe _ 2 (l)- 

Therefore, the decoder input and its output soft symbol estimates, x Ck and Xb k , for the code, 
generated by gFBc(x) = x 2+ x+1 , have the following relationship 

%b k — x Ck x Ck _ 1 x Ck _ 2 (5) 

By taking the natural logarithm of both sides of the above equation, we get 

lnx bk = lnx Ck + /nx Cfc _ 1 + lnx Ck _ 2 (6) 

We define the decoder with the input and output being the logarithm of the soft symbol 
estimates (SSE) of the coded symbols and SSEs of the information symbols, as the Log-domain 
soft-input- soft-output (SISO) decoder. As shown in Fig. [2l the SISO decoder can be implemented 
by adding a logarithm module and an exponential module at the front and rear end of the log- 
domain SISO decoder, respectively. 
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Fig. 2. The relationship of a SISO decoder and its Log-domain SISO decoder 
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Fig. 3. The Log-domain SISO forward decoder implemented by using its dual convolutional encoder 



Based on Eq. [6l log-domain SISO forward decoding of the code gFBc(x) = x2 + x+1 can 
be implemented by using the convolutional encoder, generated by the generator polynomial 
1/ ' 9fbc{%) = + x + 1, as shown in Fig. [3] Here the addition operation in the encoder is not 
carried out in the binary domain as in conventional convolutional encoders, but in the complex 
number domain. 

Eq. [6] and Fig. [3] reveal a duality relationship of the binary encoder and SISO forward decoder 
of a rate-1 feedback only convolutional code. This can be generalized to any FBC codes as 
summarized in the following theorem. 

Theorem 1 - Forward decoding duality of a feedback only convolutional (FBC) code: 

For a FBC code, generated by a generator polynomial gFBc(x) = l/q(x), we define its dual 
encoder as the encoder with the inverse generator polynomial of gFBc(x), given by qFBc(%) — 
1/ ' 9fbc{ x ) — l( x )- Then the log-domain SISO forward decoding of the FBC code can be simply 
implemented by its dual encoder in a complex field. This duality property is shown in Fig. 0] 
Proof: See Appendix A. 
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Fig. 4. Duality of a FBC encoder and its Log-domain SISO forward decoder 



Si 






S 2 






Fig. 5. The encoder and trellis of gFFc(x) = x 2 + x + 1 



B. Forward decoding of feed-forward only convolutional (FFC) code 

In this sub-section, we investigate the forward decoding of a FFC code. As will be shown in 
the following example, the duality property shown in Theorem 1 does not apply to such codes. 

Example 2: We consider a FFC code with the generator polynomial of gFFc(x) = x 2 + x + l 
for which the trellis diagram and encoder are shown in Fig. |5j 

Let lnx bk represent the output of the log-domain dual encoder, generated based on Theorem 1, 
with the generator polynomial of qFFc( x ) = l/gFFc( x ) = 1/ (x 2 + x + 1). The Table I compares 
lnxb k with the actual forward MAP decoding soft output lnxb k ■ Their differences are highlighted 
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TABLE I 

Comparison of the dual encoder output calculated based on Theorem 1 lnx bk with the actual forward 

MAP decoding soft output lnx bk 



Log soft 
input lnx Ck 


Memory Si of 
the dual encoder 


Memory S2 of 
the dual encoder 


Log soft output of 
the dual encoder lnx bk 


Desired soft decoding 
output lnx bk 


Itzx 








I/Tlx 


1/TbX (2-^ 


Ll^lX C2 


1/TbX 





I/TLX C2 ~\~ ItXXq-^ 


1/ TLX (22 ~\~ ItIX d 


ll7,X ^2 


\ XXiX C2 ~\~ ItT/X C]_ 


Itzx 


lnx C3 + lnx C2 + 
r 1 

1 lnx ci + lnx ci 1 


I/TLX (2^ ~\~ ItXX C2 




Ltix eg ~\~Itix C2 
+ lnx ci +lnx ci 


IfXX ~\~ IiTIXq-^ 


lnx Ci +lnx C3 + lnx ci + 
\lnx C2 +lnx C2 +lnx cl +lnx cl 1 


liTLX (2^ "^JiThX (2q "H liYbX (2-^ 




lnx Ci +lnx C3 + lnx ci 
+lnx C2 +lnx C2 +lnx ci +lnx ci 


lnXc 3 +lnx C2 
+ lnx ci +lnx ci 


lnx C5 +lnx Ci + lnx C2 +lnx cl + 
1 lnx C3 +lnXc 3 +lnx C2 +lnx C2 1 

L J 

1 +lnxc 1 +lnx ci +lnx ci +lrixc 1 < 


liTLX (2c j ~\'lThX (2^ "H 1/TbX CQ^^f^X (2-^ 













in the dashed-line boxes. 

From the above table, we can see that the soft outputs of the dual encoder, generated from 
Theorem 1, lnx bk are different from the actual forward MAP decoding soft outputs lnxb k when 
k > 2. This is because the recursive structure of the dual encoder qFFc(%) an d the complex 
field addition operation of the dual encoder. It can be observed from the above table if the input 
to the dual encoder is the binary symbol and addition in the encoder is a module-2 addition, 
as in the conventional binary encoder, the difference terms shown in the dotted-line-boxes will 
become zero and the dual encoder output will be equal to the actual decoding output. However, 
the inputs to the dual encoder are the logarithms of the soft inputs, which are complex numbers, 
and the addition in the dual encoder is done in the complex-number domain, which causes the 
differences between lnx bk and lnx bk . We can observe from the table that the difference terms 
come from the common terms of the shift-register contents Si and S 2 in the dual encoder. If we 
can change structure of the dual encoder by multiplying both the numerator and denominator by 
a common polynomial, without changing its actual generator polynomial, such that the encoder 
contents do not share any common elements at any time instant, then the difference between 
lnx bk and lnx bk will disappear and the dual encoder output will be equal to the actual MAP 
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Fig. 6. The modified dual encoder of g FF c{x) = x 2 + x + 1 

TABLE II 

Comparison of modified dual encoder output lnx bk with the actual forward MAP decoding soft output 

lnxt k 



Log soft 


Memory Si 


Memory S2 


Memory S3 


Log soft output of 


Desired soft decoding 


input lnx Ck 








the modified dual encoder lnx), k 


output lnxt k 


lnx ci 











lnx ci 






lnx ci 








lnx C2 + lnx ci 


lnx C2 + lnxci 






lnx ci 





LtIjCc 1 Itlx 


lnx C3 + lnx C2 


lnx C4 


lnx ca 




lnx ci 


lnx Ci +lnx C3 + lnx ci 


lnx C4 +lnx C3 + lnx ci 


lnx C5 


lnx C4 +lnx cl 






lnx cs +lnx Ci + lnx C2 +lnx cl 


lnx C5 +lnx Ci + lnx C2 +lnx cl 















forward decoding output. 

In Example 2, if we multiply both the numerator and denominator of the dual encoder generator 
polynomial q(x) by (1 + x), then we have 

Viz) = / V1 , 7 = T- 3 (7) 

Fig. [6] shows the encoder with the polynomial in Eq. ©. 

The Table HI] shows the output of the modified dual decoder and actual MAP forward decoding 
output. We can see that the soft outputs of the modified dual encoder are exactly the same as 
the actual MAP forward decoding outputs. 

We can prove that for any FFC codes, we can always find a modified dual decoder to implement 
a MAP forward decoder without changing its actual generator polynomial. This is summarized 
in Theorem 2. 

Before we present the new theorem, we first define a minimum complementary polynomial. 
For a given polynomial a(x) = x n + ■ ■ ■ + a\X + 1, we define the minimum complementary 
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polynomial as the polynomial of the smallest degree, 

z(x) = x l + zi_ix l ~ Y H h z x x + 1 (8) 



such that 



a(x)z(x) = x n+l + 1 (9) 



Since a(x) = x n + ■ ■ ■ + a\x + 1 always divides x 2 " -1 + 1, the minimum complementary 
polynomial of a(x) always exists. 

Theorem 2 - Forward decoding duality of a feed-forward only convolutional (FFC) code: 

For a FFC code, generated by a generator polynomial gFFc(x) = a(x), let z(x) represent its 

minimum complementary polynomial of degree I. The log-domain SISO forward decoding of the 

FFC code can be implemented by its dual encoder with the generator polynomial of 

, v z{x) _ z{x) _ x l + zi- X x l ~ x H \-z x x + l 

qFFc{x) ~ ^xjzjxj - x^Tl ^+1 (10) 

Proof: See Appendix B. 

Theorem 2 can be easily extended to a general convolutional (GC) code as shown in the 
following corollary. 

Corollary 1 - Forward decoding duality of a general convolutional (GC) code: For a 

GC code, generated by a generator polynomial gccix) = = ^t"'!" 13 ^ 1 , let z(x) be the 
degree-/ minimum complementary polynomial of a(x). The log-domain SISO forward decoding 
of the GC code can be simply implemented by its dual encoder with the generator polynomial 
of 

, , g(x)z(x) _ g(x)z(x) _ x n+l H \-hix + l 

qGc{X} ~ a(x)z{x) ~ x n+l + 1 ~ x n+l + 1 ( } 

1 | h n+ i^x n+l - 1 + ■ ■ ■ + h lX 

x n+l + i 

where g(x)z(x) = x n+l + h n+ i^ix n+l ~ l + • ■ ■ + h\x + 1 

This relationship of a binary encoder and its dual encoder is shown in Fig. U\ Corollary 1 can 
be directly derived from Theorem 2, so we skip its proof here. 
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Fig. 7. The encoder and its dual encoder of forward decoding of a general convolutional (GC) code 



III. Duality of Encoding and Backward Decoding of Rate-1 Convolutional 

Codes 

In this section, we investigate the MAP backward decoding of rate-1 convolutional codes 
and derive its dual encoder structure. Before discussing the backward decoding, we first define 
a reverse memory-labeling of a general convolutional (GC) code. Given the encoder of a GC 
code with rational generator polynomial g(x) = = , if we change the labeling 

of the k-th shift register in the encoder from Sk to S n -k, and change their respective feed- 
forward coefficient from to a n -k, k=l,2, . . . ,n, and feedback coefficients from bk to b n ~k, 
k=l, 2, . . . , n, we will derive an encoder with a new trellis. The resulting encoder is referred to 
as the reverse memory-labeling encoder of g(x). Figs. |8(a)| and |8(b)| show the encoder and the 
reverse memory-labeling encoder of g(x). 

In a BCJR MAP backward decoding, the received signals are decoded backward in a time- 
reverse order. That is, given the received signal sequence y = (yi, y 2 , ■ ■ ■ , yn), the order of signals 
to be decoded is from yx, yx-i, till y%. In order to decode the received signals backward, the 
decoder has to follow the trellis in a reverse direction. Figs. |9(a)] and |9(b)] show the encoder and 
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(a) The encoder of g(a;) = a(x)/q(x) 
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(b) The encoder of g(x) = a(x)/q(x) with reverse memory labeling 
Fig. 8. An encoder with reverse memory labeling 
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trellis of the code with the generator polynomial g(x) = ^PdTI • Fig- |9(c)| shows the backward 
trellis. For the decoder with the backward trellis in Fig. |9(c)[ the input to the decoder is at the 



right hand side of the decoder and its output is at the left hand side, which operates in a reverse 
direction of the conventional decoder. Fig. |9(d)| shows the corresponding forward representation 
of the backward trellis, where the decoder input and output are changed to the conventional 
order. The forward representation of the backward trellis can be implemented by an encoder 
shown in Fig. |9(e)[ When we compare Figs. |9(a)] and |9(e)"| it can be easily seen that the encoder 



in Fig. |9(a)| is the encoder of code g(x) = D 2+^ +1 and that in Fig. |9(e)| is its encoder with the 
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(d) The equivalent forward representation of backward trellis 




(e) The encoder corresponds to the trellis of (d) 
Fig. 9. Trellis, backward trellis and their respective encoders for the code g(x) = D ?^ 1 +1 
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reverse memory-labeling. 

This relationship of the encoders for the forward and backward trellises can be extended to 
general rate-1 convolutional codes, as shown in the following theorem. 

Theorem 3: Given an encoder with a generator polynomial g(x) = = ^n+'''.+°^+]> the 
forward representation of its backward trellis can be implemented by its reverse memory-labeling 
encoder of the same generator polynomial g(x). 

Proof: See Appendix C. 



From Theorem 2, we know that the log-domain SISO forward decoding of a given general 
convolutional (GC) encoder with a generator polynomial g(x) = ^fl can be implemented by 
its dual encoder with the generator polynomial qcc( x ) — a(a-)z(xj > where z{x) is the degree-/ 
minimum complementary polynomial of a(x). Then according to Theorem 3, the log-domain 
SISO backward decoding of the GC code can be implemented by the reverse memory-labeling 
encoder of qcc( x )- By combining Theorems 2 and 3, we can obtain the backward decoding 
duality, which is summarized in the following Theorem. 

Theorem 4 - Backward decoding duality of a general convolutional (GC) code: We con- 
sider a general convolutional encoder with a generator polynomial of g(x) = 4^ = x l~\ ""t aiX Vi ■ 
Let z(x) be the degree-/ minimum complementary polynomial of a(x). Its log-domain SISO 
backward decoding can be implemented by its dual encoder with reverse memory-labeling and 
the generator polynomial of 

q(x)z(x) _ q{x)z(x) _ x n+l H \-hiX + l 

qcc{X} ~ a{x)z(x) ~ x n+l + 1 ~ x n+l + 1 { } 

1 K+i-ix 11 ^ 1 j Vhxx 

- x n+l + I 

This duality is shown in Fig. [KB 

From Theorem 4, we can easily derive the backward decoding duality of a feed-forward only 
convolutional (FFC) code, summarized in the following Corollary. 

Corollary 2 - Backward decoding duality of a feed-forward only convolutional (FFC) 
code: For a FFC code, generated by a generator polynomial gFFc(x) = a(x) = x n + . . .+aix+l, 
let z(x) be the degree-/ minimum complementary polynomial of a(x). Its log-domain SISO 
backward decoding can be implemented by its dual encoder with reverse memory -labeling and 
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— — ffl*- 



The encoder for the code generated by 

Sgcvv „ 

^f(.T) X H 



Dual encoder with reverse 
memory-labeling of a log- 
domain SISO backward 
decoder for the code g GC (x) , 
given by 

q{x)z(x) 




Fig. 10. The encoder and its dual encoder for backward decoding of a general convolutional code 



the generator polynomial 

( z(x) z(x) _X l + Z l _ 1 X l ~ 1 \-Z x X+l 

qFFc[x) - ^jzV) - ^TT " i^+i ( 3) 

Corollary 2 can be proved in the same way as Theorem 4, so we skip the proof here. 

For a feedback only convolutional (FBC) code, we can prove that backward decoding does 
not contribute to the MAP calculation. The BCJR MAP decoding is exactly the same as the 
forward decoding. This is summarized in the following Theorem. 

Theorem 5 - Decoding duality of a feedback only convolutional (FBC) code: For a FBC 
code, generated by a generator polynomial Qfbc{x) = l/q(x), the MAP forward decoding is 
in fact equivalent to the BCJR MAP decoding. Its log-domain SISO decoder can be simply 
implemented by the dual encoder for the MAP forward decoding with the inverse generator 
polynomial of g F Bc(x), given by q F Bc(x) = l/g FBC (x). 

Proof: See Appendix D. 
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From Theorem 5, we can see that the MAP decoder of a FBC code can be implemented by 
its dual encoder using shift registers. This significantly reduces the decoding complexity. 

IV. The Representation of Bidirectional BCJR MAP Decoding Based on 
Encoding and Forward/Backward Decoding Duality 

In the previous two sections, we have introduced the duality of channel encoding and SISO 
MAP forward/backward decoding. Based on the derived encoding-decoding duality properties, 
in this section, we represent the bidirectional BCJR MAP decoder by linearly combining outputs 
of the dual encoders for the forward and backward decoders. By comparing the bidirectional 
BCJR MAP decoding outputs with the forward and backward dual encoder outputs, combining 
coefficients are identified through computer search such that the resulting combined forward and 
backward dual encoder outputs are exactly the same as the bidirectional MAP decoding outputs. 
In this paper, we found the expressions of these combining coefficients for some commonly used 
4-state and 8-state FFC and GC codes. The expressions for higher-states codes can be obtained 
in the same way. 

Let us first call the dual encoder of the forward and backward decoding as the forward dual 
encoder and backward dual encoder, respectively. Let x b k and x b k represent the soft outputs 
of the forward and backward dual encoders. They can be calculated based on Theorems 1-4 
in Sections II and III, respectively. Then the BCJR SISO MAP decoder output, x bk , can be 
represented as the following linear combination of forward and backward dual encoder outputs, 

x bk = Uk~Xb k + ^JkX~b k , (14) 

where w£ and tj~u are the combining coefficients in real domain applied to the forward and 
backward dual encoder outputs, respectively. The combining coefficients for some 4-state and 
8-state codes are shown below. 

A. The GC code [5/7]$, generated by g G c(x) = Ji£.+i 

The combining coefficients for the 4 states GC code [5/7] 8 can be calculated as 



ujI — tr k 



l + O/^J/4, for k = 1,3,5,..., (15a) 
l + E/x 2 J/A, for k = 2,4,6, .... (15b) 
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where x Ck is the soft symbol estimate of the received encoded symbol c k , where k = 1, 2, K, 
K is the frame of codeword, and 

\K/2\ 

0= U (16) 

i=i 

\K/2\ 

E= ft x C2l . (17) 
1=1 

where \x~\ is the ceiling operation, representing the smallest integer not less than x. 



B. The FFC code [5] 8 , generated by gFFc( x ) = x 2 + 1 

The combining coefficients for the 4-state code [5] 8 can be calculated as, 

r(l + £)/4, for k = 1,3,5,..., (18a) 

U k = = S 

t(l + 0)/4, for k — 2,4,6,.... (18b) 
where the variables O and E are defined in (fl~6l) and (flTl). 



C. 77ze FFC code [7] 8 , generated by gFFc(x) = x 2 + x + 1 
Let us first define, 

1(1,1 — ^C 3! _ ( 3_ d>3 ; (19) 

where d = 0,1,2; / = 0, 1, 2, 3, \K/3] , and (Z)jv represents the modulo operation of 
Z modulo N. We will assume I d i = 1, if the index of x c in Eq. (fl~9~l) is less than 1. 
We further define the intermediate variables, P k and Q k as, 

\K/3] 

p k = n j <^-2>a^ ( 2 °) 



z=o 



LU=[A/3]+l J (fc-l>3,i 

1 lj=0 i (fc-l>3j 



J 1 / ■ 

where [x] is the round-off operation of x, representing the integer closest to x. 
Then the combining coefficients for the 4-state code [7]s can be calculated as, 



rf k = (l + P k Q k )/4, (22) 
tJ- k = (l + P k /Q k )/A. (23) 
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D. The FFC code [17] s (8 states), generated by gFFc( x ) = x 3 + x 2 + x + 1 
Let us define 

M d ,i = x C4l _ {4 _ dU , (24) 

where d — 0, 1, 2, 3 and I — 0, 1, 2, 3, [if/4] . We assume, M dj ; = 1, if the subscript of x c in 
the above equation is less than one. We also define the intermediate variables, R k , Y k and Z k 
as, 

\K/4] 

R k = Yl M {k _ 3M , (25) 

1=0 

1 li=[fc/4]+l M 0.» .... 
_ [U=[ k /A]+l M 2^ 

Then the forward and backward combining coefficients for the 8-state code [17] 8 can be 
calculated as 



fa 



[l + R k R {k)4+1 + (R k + R {k)4+1 )Y k )/8 for A; = 1,4,5,8,9,..., (28a) 

(1 + R k R {kh+1 + (R k + R {k)4+1 )Z k )/8 for k = 2, 3, 6, 7, 10, 11.., (28b) 

{l + R k R {k)4+1 + {R k + R {k)4+1 )/Y k )/8 for k — 1,4,5,8,9, .., (29a) 

(1 + R k R {k)4+l + (R k + R {k)4+1 )/Z k )/8 for A; = 2, 3, 6, 7, 10, 11.... (29b) 



£. The GC code [15/13] 8 , generated by g GC (x) = X J^ X ^ 

The forward and backward combining coefficients for the 8-state code [15/13] 8 can be calcu- 
lated as 

ut = (1 + G k G {k)7+1 W k W k+2 + G k G {k)7+3 W k W k+5 /x Ck + 

G {k)7+1 G {k)7+3 W k+2 W k+5 /x Ck )/8, (30) 



u& = (1 + G k G {k)T+1 /(W k W k+2 ) + G k G {k)7+3 /(W k W k+5 x Ck ) + 

G \ k ) 7 +iG \ k)7 +zl '(W k+2 W k+5 x Ck )) /8. (31) 
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where 

\K/T\+1 

G k = J] iV (fc _ 3 ) 7 ,,, (32) 

1=0 

Wfc " n r(M*+i>r)/ri Nt [ ' (33) 

(34) 

and 

N d ,i = x C7(l _ 1)+d , d = 0,1, 2, 3, 4, 5, 6, Z = 0, 1, 2, 3, \K/7] + 1 (35) 

where Njj = 1, if the subscript of x c is less than one and greater than K. 

To understand how to calculate the combining coefficients described in the above equations, 
let us now look at the following example. 

Example 3: Calculation of combining coefficients for the [7] 8 FFC code 

From (fT9b we have, 

1; ^<C3 ; 1 %Cg j ^C12 1 ^Cj{) (36) 

I\i 1, x ci , x C4 , x C7 , x ClQ , — %ck-2i (37) 

-^2,Z 1) &C2 1 1 %C$ , -^CH 5 — ^Cjf — 1 " (38) 

(39) 

From (1201) . if we define 

\K/S\ \K/S] \K/S\ 

& = n v & = n j m & = n /2 >< 

«=0 1=0 1=0 

Then can be calculated as 

Pk = 6, fo, 6, 6, 6, 6,6, 6, 6, ■-, /or k = 1, 2, 3, 4, 5, 6, 7, 8, 9, (40) 

From (|40l ), we can see that, P fc actually repeats every three terms. 
Similarly, from (1211) . we have 



iU=l J 0,i _ , n _ 1U=2 J l,» _ / rt _ Hi=2 L 



2,i 

' 3, 



_ 11«=1 u >' _ /• /O _ 1-U=2 -M* _ J- r\ — A -1-1=2 *i» _ /■ 

- -ppo — t — ~~ ^' ^ 2 ~~ TP 1 r — ~~ ^ 3 ~~ Tri r — ~~ 

lli=0 J 0,i llj=0 J l,J Ilj=0 J 2J 

V4 - -j=pi - U/ J 0,1) V5 - ~=2 ; _ S2/-'l,2, V6 - ~j=f2 " ~ W J 2,2, — 

llj=0 J 0J llj=0 J lJ llj=0 J 2,i 



21 



Thus, 

Q k = Ci, c 2 , c 3 , Ci// 2 i, c 2 /a 2 2 , c 3 // 2 2 2 , Ci/(ioA 2 i G/iiyis), ■- (4i) 

That is, Qk also follows a pattern over time. To calculate the Qk for k > 3, we just need to do 
small modifications of Qks- 

V. Simulation Results 

In this section, we provide the simulation results. All simulations are performed for the BPSK 
modulation and a frame size of K=\2S symbols over AWGN channels. 

Figs. CEUto EED show the bit error rate (BER) performance of various 4-state and 8-state GC 
and FFC codes, where the curve 'Forward dual encoder + Backward dual encoder' refers to the 
direct summation of forward and backward dual encoder outputs, i.e., x b k +x & fc , and the curve 
'Linear combination of forward and backward dual encoder' refers to the optimal combining 
described in Section IV, i.e., uJk x b k +^k £ b k - 

From figures, we can see that direct summation of forward and backward dual encoder outputs 
has about 0.9dB performance loss when compared to the bidirectional BCJR MAP decoding 
for the GC code [5/7] 8 at the BER of 1CT 5 . This performance loss is reduced to around 0.2dB, 
OAdB, and 0.5dB for [7] 8 FFC, [17] g FFC, and [15] 8 FFC codes, respectively and increased to 
about l.ldB for the [15/13] 8 GC code. However, when we apply the linear combination detailed 
in Section IV to the forward and backward dual encoder outputs, their performance is exactly 
the same as the BCJR MAP decoding. One particular point needs to be noted is that for the FFC 
code [5] 8 the direct summation of forward and backward dual encoder outputs has the same 
performance as the MAP decoding, so no linear combination is actually required. 

VI. Conclusions 

In this paper, we revisited the BCJR SISO MAP forward and backward decoding process for 
the rate-1 convolutional codes. Dual encoder structures of forward and backward decoding for 
three different classes of rate-1 convolutional codes are derived. The input to the dual encoder 
is the logarithm of soft symbol estimates of the coded symbols obtained from the received 
signals, and the dual encoder output produces the logarithm of the soft symbol estimates of 
the information symbols. For the general convolutional (GC) codes, generated by a generator 
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Fig. 11. BER performances of [5/7]s GC code over AWGN channels 
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Fig. 12. BER performances of [5]s FFC code over AWGN channels 
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Fig. 13. BER performances of [7]s FFC code over AWGN channels 
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Fig. 14. BER performances of [15/13]s GC code over AWGN channels 
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SNR in dB 



Fig. 15. BER performances of [17]8 FFC code over AWGN channels 
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Fig. 16. BER performances of [15]s FFC code over AWGN channels 
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polynomial gGc( x ) — ^r\, the forward and backward decoding can be implemented by their 
corresponding dual encoders, which are generated by the polynomial, fjfyjjfy, where z(x) is the 
minimum complementary polynomial of a(x). The feed-forward only convolutional (FFC) code 
is just a special case of GC code, so it has the same dual encoder structures as the GC code. 
The derived duality property significantly reduced the the computational complexity of MAP 
forward and backward recursion from exponential to linear. Similarly, the bidirectional MAP 
decoder of GC and FFC codes can be implemented by linearly combining the outputs of dual 
encoders for the forward and backward decoding. For a feedback only convolutional (FBC) code 
9fbc( x ) — ^y> the bidirectional MAP SISO decoder is equivalent to the dual encoder for the 
forward decoding, with the generator polynomial q(x). 

In this paper, we have only focused on a class of rate-1 convolutional codes. Its significance 
is mainly as component codes in concatenated coding schemes, such as turbo coding. Actually 
the encoding-decoding duality derived in this paper can also be applied to other codes and other 
applications. For example, the transmission of digital signals in the presence of inter-symbol 
interference (ISI) can also be represented by a convolutional encoding process. The channel 
transfer function of a ISI channel can be represented by a rate-1 convolutional encoder. Thus the 
encoding-decoding duality properties can also be directly applied to facilitate the MAP channel 
detection in ISI channels. With the logarithm of the soft symbol estimates of ISI channel outputs 
as the input to the dual encoder, the output of dual encoder will produce the MAP detection 
output of ISI channels. Similarly, these duality properties should exist for other linear codes, 
which can be represented by a trellis diagram. 



A. Proof of Theorem 1 

Let us consider a feedback only convolutional (FBC) code, generated by a generator polyno- 
mial 



its encoder is shown in Fig. 0] Let Si(k),i = 1, . . . , n represent the state of memory i at time 
k. Then according to Fig. |4] we have 



VII. Appendix 



9fbc{x) = l/q(x) = l/(q n x n H qix + 1) 



(42) 



n 




(43) 
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S p (k) = S p - 1 (k-l),p>2, (44) 

where all summations are done in GF(2). 

We can rewrite the above equation as follows 

n 

b k = c k ®^q l S l {k-l) (45) 
i=i 

S 1 (k) = c k ,S p (k) = S p ^(k - l),p > 2 (46) 
S p (k) = c k „ p (47) 

where we assume that c k = for k < 0. 

Based on the above equation, we can derive the following binary decoder structure, where the 
input is the codeword symbol c k and the output is b k . 

Let Psjtyiw) denote the probability of memory Si(k) = w and a k (m) denote the probability 
of state m at time k. Let (mi, • • • , m n ) be the n-dimensional binary representation of m and 
(m' l5 • • • ,m' n ) be the binary representation of ml. At time k, with input c k , the state transits 
from (m' 1? • • • , m' n ) to (m 1; m 2 , • • • , m n )= (c fc , m' 1; • • • , m^) . Then we have 

n n 

®k(m)=Y[Ps l (k)(m i ) = P{c k = mi) JJP 5i(fc )(mi) (48) 

i=l i=2 



=P(c k = m l ) PsAk-i^m'jYlPs^irrii) 

m' n =0,l i=2 

= E (fl p 8 i (k-i)(m , j )\p(c k = m 1 ) 

m' n =0,l \j=l J 

~E a fc-i( m ')7fe( m/ > m ) ? 

m' 

where a fc _i(m') = [1 Ps j (k-i)(m' j ),^ k (m',m)=P(c k = mi) and mj=m 3 - + i, for j=l, 2, . . . , n-1. 
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The APP of bk = w can then be calculated as 

n 

p bk (w)=p(b k = w\y)= t[Ps i {k-i)(m' j )P(c k = m 1 ) (49) 

(m',m)eU(b(k)=w) j=l 
n 

E II^(»-i)K)P(Cfc = Tn 1 ) 

mi,mj,-,m^ mi® 22 gjm^w 
3=1 

n 

= ^ J^P(cfc_j = m'j)P(c k = mi) 

n j = l 

mi,m',,--- ,m' n , mi® 22 <lj m 'j= w 
3=1 



j=0 



3=0 

where m' = mi and qo = 1. 

Let L(6fc) represent the LLR of bk. From Eq. (|49l we can easily derive 

L(6 fc ) = L^g jCfe _^ (50) 

Following the L-sum theory [7], the right-hand side of (|50l) can be expanded as 

n 

/ n \ 1 + IJ tanh(L (^c fc _j) /2) 

L ^ = ln ^ < 51 ) 

Vi=o / 1 - T] tanh(L (qjCk-j) /2) 

j'=o 

where tanh(x/2) = ^vj. 

Then by using the following relationship between the LLR and soft symbol estimate, 

e L(b k ) _ l 

= e L( bk ) + i = tanh ( L ( b (k)m (52) 
L{b k )=\n\±^, (53) 

(|5Q|) can be further written as 

n n 

L(6 fc ) = L [Y, Wis = ^ ? = In ? , (54) 

Vi=o / 1-nWi 1 II (-'V, ,)" 

j=o j=o 

where x qjCkj denotes the soft symbol estimate of symbol qjCk-j. Obviously x qjCk _. = 1 when 

qj = and x qjCk _. = x c , k -j when qj = 1. Thus x qjCk _. = (x^)^. 
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By substituting ([54]) into (1521) . we get 

n 
j=0 



By taking the logarithm on both sides of (1551) . we have 

n 

In x bfc = ^ g-j In x Cfe _ :? (56) 

Therefore, the log-domain SISO forward decoding of the FBC code can be simply implemented 
by its dual encoder, generated by the generated polynomial qFBc(%) — ^-/9fbc{%) — QnX n + 
■ ■ - qix + l. 

This proved Theorem 1. 

B. Proof of Theorem 2 

Let us first examine the forward binary decoding. Based on the code generator polynomials, we 
can easily derive the binary decoder of codes generated by a(x) and ^fj^y^=^^^> as shown in 



Fig. |17(a)] and |17(b)| respectively. As can be seen from these figures, the binary decoder of each 
of these two codes is equivalent to the encoder generated by its respective inverse polynomial. 

Let (mi, ■ • • ,m n ) and (m[, ■ ■ ■ ,m' n ) be the n-dimensional binary representation of m and 
m! . Let (ui, ■ ■ ■ ,u n+ i) and (u' x , • ■ ■ ,u' n+l ) be the (n + /)-dimensional binary representation 
of u and u' . Assume that at time k, with input Ck, the state transits from (m[,- ■■ ,m' n ) to 



(mi,m 2 ,--- , m n ) in the binary decoder of Fig. 17(a) and transits from (u\,-- - ,u' n+ i) to 



,u n+ i) in 17(b) For a binary input sequence c = (ci,c 2 ,--- ,ck), it is well known 



that the polynomials and generate the same codeword. We thus have 

n— 1 n+l—1 

bk = aim>i + m 'n + C k= ^ Z 'J U 'i + U ' n + l + ° k ^ 

i=l j=l 

nij = m'j_i, and Uj = u'j_i,j > 2 (58) 
Then by following similar calculation in Appendix A, we have 

L(b k ) = L ^2 a i m 'i + m 'n + °k \ (59) 

tn+l-l \ 
^2 Z i U 'j + U 'n+l + °k I (60) 
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(a) The binary decoder of FFC code generated by a(x), which is equivalent to an encoder 
generated by l/a(x) 







* • * 









(b) The binary decoder of FFC code generated by 



a[x)z(x) _ + ' + 1 



z(x) 



z(x) 



, which is equivalent 



to an encoder generated by = 
Fig. 17. The binary dual encoder of a FFC code 



When the terms in the summation of the right-hand side in (1591 ) and (|60l) are statistically 
independent, we can use the L-sum theory to further expand these two equations. However, we 
can easily check that the terms m'^ i — 1, • • • , n, in (|59l , are not independent. Now let us prove 
that u'i, i = 1, • • • , n + I are statistically independent random variables. 

When < k < n + /, the state u[, % = 1, • • • , n + 1, at time k, is given by 



u' { = 0, k < i and u[ = Ck-i, k > i 
When k > n + I, the state u[, i — 1, • • • , n + 1, at time fc, is given by 

L*/(n+i)J 

^ Ck—pi 
p=0 



u' 



(61) 



(62) 



where [xj denotes the largest integer not greater than x. 

From (loTT) and (l62l) . we can see that u[, i = 1, • • • , n + 1, are statistically independent random 
variables at any time instant k. 
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Fig. 18. The SISO decoder structure, implemented with the encoder with the generator polynomial of qFFc( x ) = 



Since lij, z = 1, • • • , n + Z are statistically independent random variables, we can use the 
L-sum theory [7] to expand the right-hand side of (loTI) . By following a similar calculation as in 
Appendix A, we can obtain the following equation 

n+l-l 

ln *b k = Yl z j £ w 3 + x u > n+i + x Ck (63) 

3=1 

and, 

x Uj = x u i._^j > 2, (64) 

where Xb k , x Uj , x u >. and x Ck denotes the soft symbol estimate of symbol b k , Uj, u'p and Ck, 
respectively. Based on (1631) and (|64l) . we can derive the SISO decoder structure, shown in Fig. 
Il~8l implemented with the encoder with the generator polynomial of 

qFFc{x) = ^(x)4x) = ^TT (65) 

This proves Theorem 2. 
C. Proof of Theorem 3 



Assume that the encoder with the generator polynomial g(x) in Fig. |8(a)| transits from the 
state (m^, m' 2 , ■ ■ ■ , m' n ) at time k-1 to the state (mi, m 2 , • • • , m n ) at time k with input bj-, then 
we have 

n-1 

mi = b k + ^ g p mp + m' n , m p = m^, p > 2 (66) 
p=i 
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and the corresponding trellis output at time k is given by 

n—1 m—l 

c(k)=b k + 2J q P m' p + m' n + 2J a P ^ p + ^ m (67) 
p=i p=i 

n—1 m—l 

=b k + 1p m ' P + a P m P + m 'n + m 'm 

p=l p=l 

To prove Theorem 3, we now only need to prove that with input b k its reverse memory-labeling 
encoder transits from the state (m 1; m 2 , • • • , m n ) at time k-l to the state (m[, m' 2 , ■ ■ ■ , m' n ) at 
time k and generate the same encoder output. 

Now let us consider the reverse memory-labeling encoder with the generator polynomial g(x) 
in Fig. |8(b)| With the state (mi,m 2 , ■ ■ ■ ,m n ) at time k-l and input b k , the state at time k of 
the reverse memory -labeling encoder is given by 

n—1 n—1 n—1 

S n (k) = b k + m 1 + Y % m P+i =h + b k + Y Iv^v + m « + Y q P m 'p =m ' n ^ 
p=i p=i p=i 

S p (k) = Sp+i(k) = m p+1 = m' p (69) 

where in the step (a) of (|68l ) we have used Eq. (|66l ). 

The output of reverse memory -labeling encoder at time k is given by 

m—l m—l n—1 

c(k) = m n+1 + Y a p m p+1 + m 1 =m' m + Y a v m 'p + h + Y ^ m P + m 'n ( 70 ) 
p=i p=i p=i 

where we have used Eq. (l66l) in the last step of calculation. 

From (I68H70I) . we can see that with input b k the reverse memory-labeling encoder transits 
from the state (mi, ■ ■ ■ , m n ) at time k — 1 to the state (m' 1; m' 2 , ■ ■ ■ , m' n ) and generates the 
same encoder output as the encoder with the generator polynomial g(x). 

This proves Theorem 3. 

D. Proof of Theorem 5 

To prove Theorem 5, let us first examine the backward decoding of a FBC code. At the 
encoder of a FBC code in Fig. 5, with input b k , the state transits from (m[, ■ ■ ■ , m' n _i, m' n ) at 
time k — 1 to (m 1; m 2 , • ■ ■ , m n ) = (c^, m[, ■ ■ ■ , m' n _^\ at time k, where c k is the encoder output. 
The state transition is shown in the Fig. [191 where a, d, w=0 or 1, a = 1 — a, d = I — d, and 

w = 1 — w. 
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Fig. 19. The backward decoding trellis transition of a FBC code 

Then we apply the BCJR backward decoding as follows, 

(1) At time K, we have (3 K (0) = 1 and K (m) = for m ^ 0; 

(2) At time K-l, we have 

(5 K -i{m) = p Cjf (0) for m — 0, 1, and (3 K _i(m) — for m ^ 0, 1; 

(3) At time fT-2, we have 

&K-2(rn) = p CK _ 1 (0)p CK (0), for m = 0, 1,2,3, and f} K - 2 (m) = 0, for m^0, 1,2,3; 

(4) At time — 1> , < v < n, we have 

v-l 

p K -v(m) = Up CK _X0)Jorm = 0,l--- , T - 1, and p K -v(m) = 0, for m ^ 0, 1 • • • , T- 1; 

i=0 

At time K — t, t > n , we have 

n-l 

pK-t(m) = J] Pc K - t (°)> for a11 m - 

From the above equation, we can see that fikijn) is the same for all states when k < K — n. 
Therefore, the backward decoding does not have any contribution in the probability calculation 
of the BCJR decoding. This proves that the BCJR forward decoding is exactly the same as the 
BCJR MAP decoding for the FBC codes. 
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